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Bose-Einstein correlations of identical charged-pion pairs produced in 
hadronic Z decays are analyzed in terms of various parametrizations. The 
r-model with a one-sided Levy proper-time distribution provides a good 
description, enabling the source function to be reconstructed. 



1. Introduction 

In particle and nuclear physics intensity interferometry provides a direct 
experimental method for the determination of sizes, shapes and lifetimes of 
particle-emitting sources (for recent reviews see [H [21 [3]). In particular, 
boson interferometry provides a powerful tool for the investigation of the 
space-time structure of particle production processes, since Bose-Einstein 
correlations (BEC) of two identical bosons reflect both geometrical and 
dynamical properties of the particle radiating source. 

For our analysis we use a sample of about 500 thousand two-jet events, 
selected by the Durham algorithm [6] with y cut = 0.006, from e + e~annihi- 
lation data collected by L3 at a center-of-mass energy of 91.2 GeV. 



2. Parametrizations of BEC 

The two-particle correlation function is defined as: 



t> i \ P2(Pl,P2) m 

R2{Pl,P2) = f w r , (1) 
Pl{Pl)Pl{P2) 
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where P2(pi,P2) is the two-particle invariant momentum distribution, p\(pi) 
the single-particle invariant momentum distributions and pi the four-mo- 
mentum of particle i. Since we are only interested in BEC, the product 
of single particle densities is replaced by the so-called reference sample, 
/?o(pi,P2), the two-particle density that would occur in the absence of Bose- 
Einstein interference. Here we use mixed events as a reference sample [7j. 

After some assumptions [U [2] , this two-particle correlation function is 
related to the Fourier transformed source distribution. In this case 

R2(P1,P2) = 1 + \HQ)\ 2 , (2) 

where Q is the invariant four-momentum difference, Q = \J — (pi — P2) 2 and 
f(Q) is the Fourier transform of the density of the source, f(x). 



2.1. Gaussian distributed source 
The simplest assumption is that the source has a symmetric Gaussian 
distribution, in which case }{Q) = exp (ipQ — i^Lj an d 

i? 2 (Q) = 7 [l + A exp (-(RQ) 2 )} (1 + 8Q) , (3) 



where the parameter 7 is a constant of normalization, A is an intercept 
or incoherence factor, which measures the strength of the correlation, and 
(1 + 5Q) is introduced to parametrize possible long-range correlations not 
adequately accounted for in the reference sample. 

A fit of Eq.([3|) to the data results in an unacceptably low confidence 
level [7] from which we conclude that the shape of the source deviates from 
a Gaussian. The fit is particularly bad at low Q values. 



2.2. Levy distributed source 

Adopting Nolan's S(a,f3 = 0, 7, 5;1) convention [5] for the symmetric 
Levy stable distribution with rescaling of the scale parameter 7 to R and the 
location parameter 5 to xq, the Fourier transform (characteristic function) 
f(Q) has the following general form: 

f(Q) = ex V (iQx -\RQ\ a ). (4) 

The index of stability, a, satisfies the inequality < a < 2. The case a = 2 
corresponds to a Gaussian source distribution. For more details see [5]. 
Then R2 has the following, relatively simple form [9]: 

R 2 (Q) =7[l + Xexp(-(RQ) a )](l + 5Q). (5) 

After fitting Eq.([5]) to the data it is clear that the correlation function is far 
from Gaussian: a ~ 1.3. The confidence level, although improved compared 
to the fit of Eq. ([3|) , is still unacceptably low [7] . 
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3. The r-model 

A model of strongly correlated phase-space was developed [12] to ex- 
plain the experimentally found invariant relative momentum Q dependence 
of Bose-Einstein correlations in e + e~ reactions. This model also predicts 
a specific transverse mass dependence of R 2 , that we subject to an experi- 
mental test here. In this model, it is assumed that the average production 
point x^ 1 of particles with a given momentum k^ is given by 

= dk". (6) 



In the case of two-jet events, d = where r = \Jt 2 — k\ is the Ion 



gitudinal proper-time and m t = ym 2 +p 2 is the transverse mass. The 
second assumption is that the distribution of x tJ, (k tM ) about its average, 
5^(x IJ '(k^) — x^(k^)), is narrower than the proper-time distribution. Then 
the emission function of the r-model is 

S{x,k) = dTH(T)S A {x-dk)Nx(k), (7) 
Jo 

where H(t) is the longitudinal proper-time distribution, the factor 5a(x — 
dk) describes the strength of the correlations between coordinate space 
and momentum space variables and N\{k) is the experimentaly measurable 
single-particle spectrum. In the plane-wave approximation the Yano-Koonin 
formula [13] gives the following two-pion multiplicity distribution: 

P2(ki,k 2 ) = J d 4 xid 4 x 2 S'(xi, /ci)5(x 2 , k 2 ) (1 + cos [{ki - k 2 ){xi - x 2 )]) . 

(8) 

Approximating the 5/\ function by a Dirac delta function, the argument of 
the cosine becomes 

(h - k 2 )( Xl - x 2 ) = -0.5(di + d 2 )Q 2 . (9) 

Then the two-particle Bose-Einstein correlation function is obtained as 

^i^Hl + ARe^jQ, (10) 

where H{oo) = J drH (t) exp^icor) is the Fourier transform of H{t). Thus 
an invariant relative momentum dependent BEC appears. 

Guided by the result of the previous section, we use a one-sided Levy 
distribution for the longitudinal proper-time density. The corresponding 
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BEC function has an analytic form [9l I10| ITT]: 

R2(Q 2 ,m t ) = 



1 + A cos h A 



m t 



m t 



exp 




where the parameter tq is the proper-time of the onset of particle production, 
At is a measure of the width of the proper-time distribution, A = tan (^p) 
and B = (1 + SQ). 

Assuming that particle production starts immediately and defining an 
effective radius, R, [TTJ R2 simplifies to 



R2{Q)=i 



1 + A cos 



(R a Q) 2a exp(-(RQy a ) {1 + SQ), (12) 



2n 



where R a is related to R by R^ a = tan(a7r/2)i? 



l 2a . 

The fit of Eq. (fl~2|) to the data is statistically acceptable [7J. The data 
are well described by the fit. For Q between 0.5 GeV and 1.5 GeV the data 
points go below the level of the long-range correlations extrapolated to lower 
Q values. These data points indicate an anti-correlation in the Q ~ 1 GeV 
region. This property of the data is well reproduced by the fitted curve, 
which also goes below unity as a result of the cosine term in Eq. (|12p , which 
comes from the asymmetric Levy assumption. 

After fitting Eq. ([TT]) for various mt intervals we find that the quality of 
the fits is statistically acceptable and the fitted values of the model param- 
eters are stable and within errors the same in all investigated mt intervals, 
confirming the m t -dependence predicted by the r-model. The r-model with 
a one-sided Levy proper-time distribution describes the data with parame- 
ters To = fm, a ~ 0.43 ± 0.03 and At ~ 1.8 ± 0.4 fm (the difference in mt 
of the two pions is required to be less than 0.2 GeV). 



4. Reconstruction of the emission function 

In order to reconstruct the space-time picture of the emitting process we 
assume that the emission function can be factorized in the following way: 

S(r,z,t)=I(r)G( V )H(T), (13) 

where J(r) is the single-particle transverse distribution, G(rj) is the space- 
time rapidity distribution of particle production, which approximately coin- 
cides with the single-particle rapidity distribution, and H{t) is the observed 
proper-time distribution. 

With these assumptions one can reconstruct the longitudinal part of the 
emission function integrated over the transverse distribution. It is plotted 
as a function of t and z in Fig. [TJ It exhibits the typical boomerang shape 



Fig. 1. Two views of the longitudinal part of the source function normalized to the 
average number of pions per event. 



with a maximum at low t and z but with tails reaching out to very large t 
and z values. 

The transverse profile, which follows from Eq. (|7|), is given by 



This equation describes the particle production in coordinate space as a 
function of the proper-time r. It describes the expansion of the source as the 
proper-time increases. The particle production probability is proportional 
to the proper-time distribution H(t). Fig. [2] shows the transverse part of 
the emission function for various proper-times. Particle production starts 
immediately, increases rapidly and decreases slowly. A ring-like structure, 
similar to the expanding, ring-like wave created by a pebble in a pond, is 
reconstructed from L3 data, as shown in Fig. 3. An animated gif file that 
shows this effect is available from |16j . 
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